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We study the classical non-relativistic two-dimensional one-component plasma at 
Coulomb coupling T = 2 on the Riemannian surface known as Flamm's paraboloid 
which is obtained from the spatial part of the Schwarzschild metric. At this special 
value of the coupling constant, the statistical mechanics of the system are exactly 
solvable analytically. The Helmholtz free energy asymptotic expansion for the large 
system has been found. The density of the plasma, in the thermodynamic limit, has 
been carefully studied in various situations. 
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I. INTRODUCTION 

The system under consideration is a classical (non quantum) two-dimensional one- 
component plasma: a system composed of one species of charged particles living in a 
two-dimensional surface, immersed in a neutralizing background, and interacting with the 
Coulomb potential. The one-component classical Coulomb plasma is exactly solvable in one 
dimension In two dimensions, in their 1981 work, B. Jancovici and A. Alastuey 

aa 

showed how the partition function and n-body correlation functions of the two-dimensional 
one-component classical Coulomb plasma (2dOCP) on a plane can be calculated exactly 
analytically at the special value of the coupling constant T = (3q 2 = 2, where f3 is the inverse 
temperature and q the charge carried by the particles. This has been a very important re- 
sult in statistical physics since there are very few analytically solvable models of continuous 
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fluids in dimensions greater than one. 

Since then, a growing interest in two-dimensional plasmas has lead to study this system 



on various flat geometries 
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5|, |6j and two-dimensional curved surfaces: the cylinder 
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the sphere jl], [lO, III, ll2l, ll3| and the pseudosphere [l_4|, ll5|, ll6|. These surface have constant 
curvature and the plasma there is homogeneous. Therefore, it is interesting to study a case 
where the surface does not have a constant curvature. 

In this work we study the 2dOCP on the Riemannian surface S known as the Flamm's 
paraboloid, which is obtained from the spatial part of the Schwarzschild metric. The 
Schwarzschild geometry in general relativity is a vacuum solution to the Einstein field equa- 
tion which is spherically symmetric and in a two dimensional world its spatial part has the 
form 

rfs 2 = (l-—) 1 dr 2 + r 2 dif 2 . (1.1) 



In general relativity, M (in appropriate units) is the mass of the source of the gravitational 
field. This surface has a hole of radius 2M and as the hole shrinks to a point (limit M — > 0) 
the surface becomes flat. It is worthwhile to stress that, while the Flamm's paraboloid 
considered here naturally arises in general relativity, we will study the classical (i.e. non 
quantum) statistical mechanics of the plasma obeying non-relativistic dynamics. Recent 
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developments for a statistical physics theory in special relativity have been made in 
To the best of our knowledge no attempts have been made to develop a statistical mechanics 
in the framework of general relativity. 

The "Schwarzschild wormhole" provides a path from the upper "universe" to the lower 
one. We will study the 2dOCP on a single universe, on the whole surface, and on a single 
universe with the "horizon" (the region r = 2M) grounded. 

Since the curvature of the surface is not a constant but varies from point to point, the 
plasma will not be uniform even in the thermodynamic limit. 

We will show how the Coulomb potential between two unit charges on this surface is given 
by — ln(|^i— z-i\j \J\z\Z-i\) where Z\ = (^/ri+V r i ~ 2M) 2 e llpi . This simple form will allow us to 
determine analytically the partition function and the n-body correlation functions at T = 2 
by extending the original method of Jancovici and Alastuey 

aa 

. We will also compute the 

thermodynamic limit of the free energy of the system, and its finite-size corrections. These 
finite-size corrections to the free energy will contain the signature that Coulomb systems 
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can be seen as critical systems in the sense explained in 




The work is organized as follows: in section [TTJ we describe the one-component plasma 
model and the Flamm's paraboloid, i.e. the Riemannian surface S where the plasma is 
embedded. In section IHIl we find the Coulomb pair potential on the surface S and the 
particle-background potential. We found it convenient to split this task into three cases. 
We first solve Poisson equation on just the upper half of the surface S. We then find the 
solution on the whole surface and at last we determine the solution in the grounded horizon 
case. In section IIV} we determine the exact analytical expression for the partition function 
and density at T = 2 for the 2dOCP on just one half of the surface, on the whole surface, 
and on the surface with the horizon grounded. In section [V] we outline the conclusions. 



II. THE MODEL 



A one-component plasma is a system of N pointwise particles of charge q and density n 
immersed in a neutralizing background described by a static uniform charge distribution of 
charge density pb = —qn b . 

In this work, we want to study a two-dimensional one-component plasma (2dOCP) on a 
Riemannian surface S with the following metric 

(2M\ 1 
1 J dr 2 + r 2 d V 2 . (2.1) 

or g„. = 1/(1 - 2M/r),g w = r 2 , and g rip = 0. 

This is an embeddable surface in the three-dimensional Euclidean space with cylindrical 
coordinates (r, <p, Z) with ds 2 = dZ 2 + dr 2 + r 2 d(p 2 , whose equation is 

Z(r) = ±2y/2M(r - 2M) . (2.2) 

This surface is illustrated in Fig. [TJ It has a hole of radius 2M. We will from now on call 
the r = 2M region of the surface its "horizon" . 



1. The Flamm's paraboloid S 



The surface S whose local geometry is fixed by the metric (11.11) is known as the Flamm's 
paraboloid. It is composed by two identical "universes": S + the one at Z > 0, and <S_ 



FIG. 1: The Riemannian surface S: the Flamm's paraboloid. 



the one at Z < 0. These are both multiply connected surfaces with the "Schwarzschild 
wormhole" providing the path from one to the other. 

The system of coordinates (r, ip) with the metric (11.11) has the disadvantage that it requires 
two charts to cover the whole surface S. It can be more convenient to use the variable 

Z 



u = = ±W 1 (2.3) 

AM V 2M K J 

instead of r. Replacing r as a function of Z using equation (12.21) gives the following metric 
when using the system of coordinates (u,<p), 



ds 2 = AM 2 (l + u 2 ) [Adu 2 + {l + u 2 ) dcp 2 ] . 

The region u > corresponds to S + and the region u < to <S_. 
Let us consider that the OCP is confined in a "disk" defined as 



(2.4) 



Q+ = { q = (r, ip) e S + \0 < ip < 2tt, 2M < r < R} . 



The area of this disk is given by 
Ar = / dS = tc 



^R(R-2M)(3M + R) + 6M 2 In 



V2M 



(2.5) 



(2.6) 
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where dS = ^fgdr dip and g = det(g^ u ). The perimeter is Cr = 2itR. 

The Riemann tensor in a two dimensional space has only 2 2 (2 2 — 1)/ 12 = 1 independent 
component. In our case the characteristic component is 

M 

R wrw • (2.T) 

r 

The scalar curvature is then given by the following indexes contractions 

TZ = R% = W\ v = 2iT% = 2g^R r wv = -™ , (2.8) 

and the (intrinsic) Gaussian curvature is K = 1Z/2 = —M/r 3 . The (extrinsic) mean curva- 
ture of the manifold turns out to be H = — ^ M/8r 3 . 
The Euler characteristic of the disk is given by 

X = _L I / RdS+ I kdl) , (2.9) 



2tt 



9Q R 



where k is the geodesic curvature of the boundary dfl~^. The Euler characteristic turns out 
to be zero, in agreement with the Gauss-Bonnet theorem \ = 2 — 2h — b where h = is the 
number of handles and 6 = 2 the number of boundaries. 

We can also consider the case where the system is confined in a "double" disk 

n R = n+un R , (2.io) 

with = {q = (r, tp) e 5_|0 < ip < 2ir, 2M < r < R}, the disk image of on the lower 
universe <S_ portion of S. The Euler characteristic of Qr is also x = 0. 



2. A useful system of coordinates 



The Laplacian for a function / is 

1 d 



A/ 



Vgg 



d_ 

dq h 



f 



] _ 2M\ d 2 1 ^ 

r J dr 2 r 2 dip 2 



1 M\ d 
- I dr- 



/ 



(2.11) 



where q = (r, ip). In appendix [A] we show how, finding the Green function of the Laplacian, 
naturally leads to consider the system of coordinates (x,<p), with 



x 



(Vu 2 + l + uf 



[2.12) 



6 



The range for the variable x is ]0, +00 [. The lower paraboloid <S_ corresponds to the region 
< x < 1 and the upper one S + to the region x > 1. A point in the upper paraboloid with 
coordinate (x, ip) has a mirror image by reflection (u — > —u) in the lower paraboloid, with 
coordinates (l/x,(p), since if 



x = (Vw 2 + 1 + uf (2.13) 

then 

- = (Vm 2 + 1 - u) 2 . (2.14) 

x 

In the upper paraboloid S + , the new coordinate x can be expressed in terms of the original 
one, r, as 



U/F+V^W (215) 

2M v ; 
Using this system of coordinates, the metric takes the form of a flat metric multiplied by 
a conformal factor 

cis 2 = ^ ^1 + ^ (rfx 2 + x 2 V) • (2.16) 
The Laplacian also takes a simple form 

A/ = j Aflat/ (2.17) 

M 2 (l + i) 4 1 ^ 

where 

A f-^l + -^L + L^L (2 18) 
(9x 2 x dx x 2 <9(/9 2 

is the Laplacian of the flat Euclidean space R 2 . The determinant of the metric is now given 
by g = [M 2 x(l +x- r ) A /4] 2 . 

With this system of coordinates (x, ip), the area of a "disk" of radius R [in the original 
system (r, ip)] is given by 

A R = —-p{x m ) (2.19) 



with 



1 



p(x) = x 1 + 8x + 12 In x (2.20) 

XX 2 



and x rn = (VR+ ^/R-2M) 2 /(2M). 



7 



III. COULOMB POTENTIAL 

A. Coulomb potential created by a point charge 

The Coulomb potential G(x, ip; xq, <po) created at (x,ip) by a unit charge at (x ,y?o) is 
given by the Green function of the Laplacian 

AG(ar, ip; x , v?o) = -2tt<5 (2) (x, <p\ x , p ) (3.1) 
with appropriate boundary conditions. The Dirac distribution is given by 

5 {2 \x, ip; x , ipo) = M 2 x( ^ +x -iy 6 ( x ~ x o)% ~ <Po) ( 3 ' 2 ) 

Notice that using the system of coordinates (x, (p) the Laplacian Green function equation 
takes the simple form 

A flat G(x, p; x , ip ) = -2tt- 5{x - x )5(<p - ip Q ) (3.3) 

x 

which is formally the same Laplacian Green function equation for flat space. 

We shall consider three different situations: when the particles can be in the whole surface 
S, or when the particles are confined to the upper paraboloid universe <S + , confined by a 
hard wall or by a grounded perfect conductor. 

1. Coulomb potential G ws when the particles live in the whole surface S 

To complement the Laplacian Green function equation (13.11) . we impose the usual bound- 
ary condition that the electric field — VG vanishes at infinity (x — > oo or x — > 0). Also, 
we require the usual interchange symmetry G(x, <p; xo, <po) = G(xo, ipo; x, ip) to be satisfied. 
Additionally, due to the symmetry between each universe S + and iS_, we require that the 
Green function satisfies the symmetry relation 

G ws (x, ^ x , p ) = G ws (l/x, p; l/x , <p ) (3.4) 

The Laplacian Green function equation (13. ip can be solved, as usual, by using the de- 
composition as a Fourier series. Since equation (13. ip reduces to the flat Laplacian Green 
function equation (13. 3p . the solution is the standard one 

°° 1 fx \ 2n 

G(x, ip; x , ipo) = ^2 ~ ( ~ ) cos [n(ip - p> )]+ 9o(x,x ) (3.5) 
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where x > = max(x, Xq) and x< = min(x, Xq). The Fourier coefficient for n = 0, has the form 

{dp In x + &o 5 x > x o 
(3.6) 
do In x + 6q , x < xo • 

The coefficients a^&o are determined by the boundary conditions that go should be con- 
tinuous at its derivative discontinuous d x g \ x=x + — d x g \ x=x - = —l/x , and the 
boundary condition at infinity V<7oU-kx> = and V<7oU->o = 0. Unfortunately, the boundary 
condition at infinity is trivially satisfied for g Q , therefore go cannot be determined only with 
this condition. In flat space, this is the reason why the Coulomb potential can have an ar- 
bitrary additive constant added to it. However, in our present case, we have the additional 
symmetry relation ( 13. which should be satisfied. This fixes the Coulomb potential up to 
an additive constant b Q . We find 

g {x,x ) = --In— + b , (3.7) 

2* X <- 

and summing explicitly the Fourier series (13. 5ft . we obtain 

G ws (x, <p; x , <p ) = - In + b ° > ( 3 - 8 ) 

V\ zz o\ 

where we defined z = xe ltp and zq = xoe t(po . Notice that this potential does not reduce 
exactly to the flat one when M = 0. This is due to the fact that the whole surface S in the 
limit M — > is not exactly a flat plane M. 2 , but rather it is two flat planes connected by a 
hole at the origin, this hole modifies the Coulomb potential. 



2. Coulomb potential G when the particles live in the half surface S + confined by hard walls 

We consider now the case when the particles are restricted to live in the half surface S + , 
x > 1, and they are confined by a hard wall located at the "horizon" x — 1. The region 
x < 1 (<S_) is empty and has the same dielectric constant as the upper region occupied 
by the particles. Since there are no image charges, the Coulomb potential is the same G ws 
as above. However, we would like to consider here a new model with a slightly different 
interaction potential between the particles. Since we are dealing only with half surface, we 
can relax the symmetry condition (13.41) . Instead, we would like to consider a model where 
the interaction potential reduces to the flat Coulomb potential in the limit M — > 0. The 
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solution of the Laplacian Green function equation is given in Fourier series by equation (13.51) . 
The zeroth order Fourier component go can be determined by the requirement that, in the 
limit M — > 0, the solution reduces to the flat Coulomb potential 



G flat (r,r') = -ln^— ^ (3.9) 

LJ 

where L is an arbitrary constant length. Recalling that x ~ 2r/M, when M — > 0, we find 

M 

go(x,x ) = -lnx> - In— (3.10) 



and 



M 

G s (x, <p; x , (fo) = - In \z - z \ - In — . (3.11) 



5. Coulomb potential G gh «//ien £/ie particles live in the half surface 5+ confined by a grounded 

perfect conductor 

Let us consider now that the particles are confined to S + by a grounded perfect conductor 
at x = 1 which imposes Dirichlet boundary condition to the electric potential. The Coulomb 
potential can easily be found from the Coulomb potential G ws (13.81) using the method of 
images 

I y — y n I I y — r_ ^ I y — r 

(3.12) 



G gh (x, if; x , ip ) = - In ^j=^L + In — %= = - In 

y/\zzo\ V\ zz o I 



z 



1 - zz 

where the bar over a complex number indicates its complex conjugate. We will call this the 
grounded horizon Green function. Notice how its shape is the same of the Coulomb potential 
on the pseudosphere 15] or in a flat disk confined by perfect conductor boundaries [4|. 

This potential can also be found using the Fourier decomposition. Since it will be useful 
in the following, we note that the zeroth order Fourier component of G gh is 

g (x,x ) = lnx< . (3.13) 



B. The background 

The Coulomb potential generated by the background, with a constant surface charge 
density pb satisfies the Poisson equation 

Av b = -27i Pb . (3.14) 
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Assuming that the system occupies an area A.r, the background density can be written 
as p b = —qN b /AR = —qn b , where we have defined here n b = N b /An the number density 
associated to the background. For a neutral system N b = N. The Coulomb potential of 
the background can be obtained by solving Poisson equation with the appropriate boundary 
conditions for each case. Also, it can be obtained from the Green function computed in the 
previous section 

v b (x, (f) = J G(x, cp; x\ cp')p b dS' . (3.15) 

This integral can be performed easily by using the Fourier series decomposition (13. 5j) of the 
Green function G. Recalling that dS = |M 2 x(l + x^ 1 ) 4 dx dip, after the angular integration 
is done, only the zeroth order term in the Fourier series survives 

. np b M 2 [ Xm . ,, / l\ 4 , . 
v b (x,ip) = — - — J g (x,x)x\l + -\ dx. (3.16) 

The previous expression is for the half surface case and the grounded horizon case. For the 
whole surface case, the lower limit of integration should be replaced by l/x m , or, equivalently, 
the integral multiplied by a factor 2. 

Using the explicit expressions for g , (13.71) . (I3.10p . and (13.131) for each case, we find, for 
the whole surface, 

v b s (x, if) = - npb ^ [h( x ) - h(x m ) + 2p(x m ) lnx m - Ab p(x m )} (3.17) 

o 

where p(x) was defined in equation (12.201) . and 



h(x) = x 2 + 16x + — + ^ + 12(lnx) 2 - 34 . (3.18) 

x x l 



Notice the following properties satisfied by the functions p and h 

p( x ) = -p(l/x) , h(x) = h(l/x) (3.19) 

and 

x x4 



p{x) = xh'(x)/2 , p'(x) =2x\l + -J (3.20) 

where the prime stands for the derivative. 

The background potential for the half surface case, with the pair potential — ln(|z — 
z'\M/2L) is 



w *PbM 2 



v?(x,tp) 



x m M 

h(x) - h(x m ) + 2p(x m ) In -^j- 



(3.21) 
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Also, the background potential in the half surface case, but with the pair potential — ln(|z- 
z'l/y^D + bo is 



irpbM 2 



h(x) 



h(x 



Finally, for the grounded horizon case, 

h np b M 2 
% {x, <P) = — 3 — 



h p(x m ) In 2b 

2 \ x 



[h(x) — 2p{x m ) lnx] . 



(3.22) 



(3.23) 



IV. PARTITION FUNCTION AND DENSITY AT V = 2 

We will now show how at the special value of the coupling constant Y = f3q 2 = 2 the 
partition function and n-body correlation functions can be calculated exactly. 

In the following we will distinguish four cases labeled by A: A = hs, the plasma on the 
half surface (choosing G hs as the pair Coulomb potential); A = ws, the plasma on the whole 
surface (choosing G ws as the pair Coulomb potential); A = hs, the plasma on the half surface 
but with the Coulomb potential G ws of the whole surface case; and A = gh, the plasma on 
the half surface with the grounded horizon (choosing G gh as the pair Coulomb potential). 

The total potential energy of the plasma is, in each case 

V A = V$ + q Vb(?i) + <?' ° A ( Xi > ^ ^' ' 
i i<j 

where (xj, ipi) is the position of charge i on the surface, and 

< = \j Pbrf(x,<p)dS (4.2) 

is the self energy of the background in each of the four mentioned cases. In the grounded 
case A = gh, one should add to V sh in (14. ip the self energy that each particle has due to 
the polarization it creates on the grounded conductor. 



A. The 2dOCP on half surface with potential - In \z - z'\ - lnM/(2L) 

1. Partition function 



For this case, we work in the canonical ensemble with N particles and the background 
neutralizes the charges: N b = N, and n = N/Ar = n b . The potential energy of the system 
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takes the explicit form 



a 2 N a 2 M a 2 

V hs = -q 2 ^\zi-z j \ + ^aJ2Hxi) + \Nln—-'LNah(x m ) 

l<i<j<N i=l 

a 2 a 2 f Xm 

+^-N 2 \nx m -^-a 2 h(x)p'{x)dx (4.3) 
2 4 J i 

where we have used the fact that dS = ttM 2 x(1 + x^ 1 ) 4 dx/2 = ixM 2 p'(x) dx/4, and we have 
defined 

™ b M 2 

a = . (4.4) 

Integrating by parts the last term of (14. 3[) and using (13.201) . we find 

a 2 N a 2 M n 2 

V hs = -q 2 ^\z l -z J \ + \aY,Kxi) + \N\n— + q -N 2 \nx m 



2 ^— ' 2 2L 2 

l<i<j<N i=l 

+ q*2 f Xm WW dx _ t Nah{Xm) . (4 .5) 
z J i x 2 

When f3q 2 = 2, the canonical partition function can be written as 

^ hs = ^^o hs exp(-/?F hs ) (4.6) 

with 

M f Xm \n(r)] 2 

- (3F^ S = -N In — - N 2 In x m - a 2 / dx + Nah(x m ) (4.7) 
2L J 1 x 



and 

N 

n?hs 



' i=l l<i<i<AT 

where A = \j2irf3h 2 /m is the de Broglie thermal wavelength. Z can be computed using 
the original method for the OCP^hi flat space , which was originally introduced in 



the context of random matrices [19|, |20|]. By expanding the Vandermonde determinant 
n«<j( z i — z j) an d performing the integration over the angles, the partition function can be 
written as 

N-l 

4 s = n m*o , (4.9) 

k=0 

where 

B N {k) = J x 2k e~ ah(x) dS (4.10) 



— / x 2k e~ ahix) p'(x)dx. (4.11) 
nb Ji 
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In the flat limit M — > 0, we have x ~ 2r/M, with r the radial coordinate of the flat space 
IR 2 , and h(x) ~ p(x) ~ x 2 . Then, Bn reduces to 

BvW-J^-jik + ^N) (4.12) 

where 7(fc + l,N) = t fc e~* dt is the incomplete Gamma function. Replacing into (14.91) . 
we recover the partition function for the OCP in a flat disk of radius R 3] 

In /V 4- 

n b A 4 4 



In = - In _ + _ - — In N + £ In 7 (fc, N) . (4.13) 

fc=i 



Thermodynamic limit R —* oo, x m — > oo ; and /ixed M 



Let us consider the limit of a large system when x m = + VR~ 2M) 2 / (2M) — > oo, 
AT — ► oo, constant density n, and constant M. Therefore a is also kept constant. In 
appendix [HI we develop a uniform asymptotic expansion of £>at(A;) when N — > oo and 
fe — > oo with (iV — k)/yN = 0(1). Let us define by 

k = ap{x k ). (4.14) 

The asymptotic expansion ( 1B15I) of <Sjv(A;) can be rewritten as 

<M*0 = ^V^S) e 2fcln£fc - ah ^) [1 + erf (e*)] 



x 



(4.15) 



2p(x fc ) iV-jfe , 
e k = y ) kJ — 4.16 
XkP'(xk) V2k 



where 



is a order one parameter, and the functions £i( e fc) an d G( e fc) can be obtained from the 
calculation presented in appendix [Bj They are integrable functions for e/. G [0, oof. We will 
obtain an expansion of the free energy up to the order In N. At this order the functions £1,2 
do not contribute to the result. 
Writing down 

N 

lnZ hs = J2 lnB N(k) - In B N (N) (4.17) 



k=0 
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and using the asymptotic expansion (I4.15p . we have 



-N In + S\ s + 5 2 hs + S, hs + — In N 
V2tt 12 



-In 



\fotX m ( 1 



1 



- 2N\iax m + ah(x m ) + 0(1) 



(4.18) 



with 



k=o \ ^ k ' 

N 

5 2 hs = ^\2k\nx k -ah(x k )} 
1 + erf (e k ) 



Sf 



k=0 

N 

k=0 



(4.19) 
(4.20) 
(4.21) 



Notice that the contribution of £i(efc) is of order one, since 2~2 k £,i(tk)/Vk ~ Jo° £,i( e ) de 
0(1). Also, £*6(e*)/*~ (1/v^T) / °° &(e)de = OQ./VN). 



5*3 S gives a contribution of order yN, transforming the sum over k into an integral over 
the variable t = e k , we have 



Jo ^ 



(4.22) 



This contribution is the same as the perimeter contribution in the flat case. 

To expand and up to order 0(1), we need to use the Euler-McLaurin summation 



formula 



21 



22] 



£ f(k) = f f(y) dy+ l - [/(0) + /(JV)] + i [f (iV) - /'(0)] + • • • . (4.23) 

fc=0 ^° 



We find 



+ ( 14a + - ) ln x m + 6 (ln x r 



(4.24) 



and 



N 2 ln x m + N ln x m — aNh(x m ) + a 1 



M^L dx -%h(x m ) + ^lnx m (4.25) 
x 2 o 
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Summing all terms in IiiZq s and those from f3F^ s , we notice that all nonextensive terms 
cancel, as it should be, and we obtain 

\nZ hs = -Nf3f B + Ax m a - C R Pj haxd + (l4ot - lnx m + 0(1) (4.26) 

where 

1 7tt 2 T 2 

^ = -2 ln ^ < 427) 

is the bulk free energy of the OCP in the flat geometry 

/^ = -yf jflniif^ (4.28) 
is the perimeter contribution to the free energy ( "surface" tension) in the flat geometry near 



a plane hard wall |5], and 

C R = 2nR = 7iM^/x m p'{x m )/2 = nhlx m + 0(1) (4.29) 



is the perimeter of the boundary at x — x r 

The region x — > oo has zero curvature, therefore in the limit x m — > oo, most of the 
system occupies an almost flat region. For this reason, the extensive term (proportional to 
N) is expected to be the same as the one in flat space fs- The largest boundary of the 
system also in an almost flat region, therefore it is not surprising to see the factor 

7hard from the flat geometry appear there as well. Nevertheless, we notice an additional 
contribution 4ax m to the perimeter contribution, which comes from the curvature of the 
system. In the logarithmic correction \nx m , we notice a —(1/6) lnx m term, the same as in a 
flat disk geometry joj, but also a nonuniversal contribution due to the curvature 14alnx m . 



3. Thermodynamic limit at fixed shape: a — > oo and x m fixed 

In the previous section we studied a thermodynamic limit case where a large part of the 
space occupied by the particles becomes flat as x — > oo keeping M fixed. Another interesting 
thermodynamic limit that can be studied is the one where we keep the shape of the space 
occupied by the particles fixed. This limit corresponds to the situation M — > oo and R — > oo 
while keeping the ratio R/M fixed, and of course the number of particles N —>■ oo with the 
density n fixed. Equivalently, recalling that N = ap(x m ), in this limit x m is fixed and finite, 
and a = TrM 2 rib/4: — > oo. We shall use a as the large parameter for the expansion of the 
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free energy. In this limit, we expect the curvature effects to remain important, in particular 
the bulk free energy (proportional to a) will not be the same as in flat space. 
Using the expansion (1B18I) of <8jv(A;) for the fixed shape situation, we have 



lnZ hs 



where now 



X In ^— + X In + S* sfixed + S 2 hs ' fixed + S* sfixed + 0(1) 



N-l 



(4.30) 



^hs.fixed = \n[x k p'(x k )} 



S. 



hs, fixed 



k=0 
N-l 

—a 

k=0 

N-l 



-ihs, fixed 



E 1 * 

k=0 



^2[h(x k ) - 2p(x k )\nx k } 

a 

erf(e M ) + erf(e fe , m ) 



(4.31) 
(4.32) 
(4.33) 



with e k)Tn and e^i given in equations (1B19I) and (1B20j) . and x k is given by k = ap(x k ). Using 
the Euler-McLaurin expansion, we obtain 

" (1 + s) 4 2(x + l) 4 



ghs, fixed 



In 



x° 



S. 



hs, fixed 



X 2 In x m — aXh(x m ) + a 2 



dx + 0(l) 



(4.34) 



dx + -h{x m ) -X\nx m + 0(1X4.35) 
x 2 



For S^' ed , the relevant contributions are obtained when k is of order V~N, where e k i is of 



order one, and when X — k is of order y X, where e kjTn is of order one. In those regions, the 
sum can be changed into an integral over the variable t = e ki i or t = e k ^ m . This gives 

2 



s. 



hs, fixed 



Aira 
n b 



X 



1 + — I +4 



/?7hard + 0(1) 



(4.36) 



with 7hard given in equation (14.281) . Once again the nonextensive terms (proportional to a 2 



in 5*2 s ' xe cancel out with similar terms in F S ' xe from equation (14.71) . The final result for 
the free energy j3F hs 

\nZ hs 



In Z hs is 



a 



-p(x m )f3fB + - [h(x m ) - 2p(x m ) \nx m ] + 



Xm (1 + x) 4 . (x + l) 4 ■ 
^ — In — ax 



x" 



X, 



n b 



X ;■ 



1+— ) + I 

•Em 



/?7hard + 0(l) 



(4.37) 



where given by (I4.27p . is the bulk free energy per particle in a flat space. We notice the 
additional contribution to the bulk free energy due to the important curvature effects [second 
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and third term of the first line of equation (14.371) ] that remain present in this thermodynamic 
limit. 

The boundary terms, proportional to y/a, turn out to be very similar to those of a flat 
space near a hard wall [23( , with a contribution /?7hardC& for each boundary at x b = x m and 
at Xj, — 1 with perimeter 

C b = vrM^SS = nMx b (l + ±-) " . (4.38) 

Also, we notice the absence of In a corrections in the free energy. This is in agreement with 
the general results from Refs. 0, Q , where, using arguments from conformal field theory, it is 
argued that for two-dimensional Coulomb systems living in a surface of Euler characteristic 
X, in the limit of a large surface keeping its shape fixed, the free energy should exhibit a 
logarithmic correction (x/6) InR where R is a characteristic length of the size of the surface. 
For our curved surface studied in this section, the Euler characteristic is x — 0, therefore no 
logarithmic correction is expected. 

4- Distribution functions 
Following [2I, we can also find the k-body distribution functions 

n Whs (qi, • • • , q*) = det[/C^(q i; q^feMi,...,^ , (4.39) 
where q« = (xj, (pi) is the position of the particle i, and 

N-l k=k p -a[h(\zi\)+h(,\zj\)]/2 

C(q, Qi) = E 1 - B m • ( 44 °) 

where = x^e %iph . In particular, the one-body density is given by 

N ~ X 2k -ah(x) 

n hs (x) = Mq, q) = E » m ■ ( 4 - 41 ) 

k=o Bn ^> 

5. Internal screening 



Internal screening means that at equilibrium, a particle of the system is surrounded by 
a polarization cloud of opposite charge. It is usually expressed in terms of the simplest of 
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the multipolar sum rules 24|: the charge or electroneutrality sum rule, which for the OCP 
reduces to the relation 



n 



(2)hs. 



q 1 ,q a )d5 2 = (JV-l)nW h8 (q a 



(4.42) 



This relation is trivially satisfied because of the particular structure (I4.39j) of the correlation 
function expressed as a determinant of the kernel /C^f, and the fact that JC^ is a projector 



^3^(q 1 ,q3)^(q 3) q 2 ) = /C^q,) 



"lis 



lis/ 



(4.43) 



Indeed, 



n 



(2)hs, 



qi,q 2 ) dS 2 



[£5(qi, qi)^ s (q 2 , q 2 ) - ^( qi , q 2 )/C^ s (q 2 , ciS 2 



-lis 



lis/ 



-ha/ 



n 



(l)hs. 



lis/ 



(iV - l)n« hs ( qi ) . 



(4.44) 



6. External screening 

External screening means that, at equilibrium, an external charge introduced into the sys- 
tem is surrounded by a polarization cloud of opposite charge. When an external infinitesimal 
point charge Q is added to the system, it induces a charge density pg(q). External screening 
means that 

J p Q (q)dS=-Q . (4.45) 

Using linear response theory we can calculate pQ to first order in Q as follows. Imagine that 
the charge Q is at q. Its interaction energy with the system is H int = Q(fi(q) where 0(q) 
is the microscopic electric potential created at q by the system. Then, the induced charge 
density at q' is 

PQ(q) = -P(M)H mt ) T = -/?Q(p(q')0(q)>T , (4.46) 

where p(q') is the microscopic charge density at q', (AB) T = (AB) — (A)(B), and (. . .) 
is the thermal average. Assuming external screening (14.45)) is satisfied, one obtains the 



Carnie-Chan sum rule 



24] 



P J<M)kq»TdS'=l. (4.47) 
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Now in a uniform system starting from this sum rule one can derive the second moment 



Stillinger-Lovett sum rule 24|. This is not possible here because our system is not homo- 
geneous since the curvature is not constant throughout the surface but varies from point to 
point. If we apply the Laplacian respect to q to this expression and use Poisson equation 

A q (p(q')<Kq)>T = -27r(p(q / )/3(q))T , (4.48) 

we find 

J pi 2 )(q\q)dS' = 0, (4.49) 

where p e (q',q) = (/KqO/Kq^T is the excess pair charge density function. Eq. (14.491) is 
another way of writing the charge sum rule Eq. (I4.42p in the thermodynamic limit. 

7. Asymptotics of the density in the limit x m — > oo and a fixed, for 1 <C x <C x m 

The formula (14.411) for the one-body density, although exact, does not allow a simple 
evaluation of the density at a given point in space, as one has first to calculate B^(k) 
through an integral and then perform the sum over k. One can then try to determine the 
asymptotic behaviors of the density. 

In this section, we consider the limit x m — > oo and a fixed, and we study the density in 
the bulk of the system l<iCi m . 

In the sum (I4.4ip . the dominant terms are the ones for which k is such that x k = x, 
with x k defined in (I4.14p . Since 1 < i < x m , the dominant terms in the calculation of the 
density are obtained for values of k such that 1 <C k <C N. Therefore in the limit N — > oo, 
in the expansion (14.151) of B^ik), the argument of the error function is very large, then the 
error function can be replaced by 1. Keeping the correction l/(12/c) from (14.151) allow us to 
obtain an expansion of the density up to terms of order 0(1/ x 2 ). Replacing the sum over k 
into an integral over Xk, we have 

n hs (x) = ^= ^ e*(**)/(x fc ) (l - - V s ) dik ( 45 °) 
J-oo V I2ap(x k )j 

with 

$>(x fc ) = 2ap(x k ) In ^ - a[h(x) - h{x k )\ (4.51) 
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and 



f( x k) = \l ^ ■ (4.52) 



We proceed now to use the Laplace method to compute this integral. The function ^(x k ) 
has a maximum for x = x k , with ty(x) = and 

= - 2 -^ (4.53a) 

q?W( x ) = -— + o(l/x 2 ) (4.53b) 

x 

$( 4 )( x ) = ^ + 0(l/x 3 ). (4.53c) 
Expanding for x k close to x and for x ^> 1 up to order 1/x 2 , we have 

n hs (x) = ^= [ + °° e -V(*)(* fc -«)»/* ( f( x ) + f( x )(x k - x) + f "^ k ~ X ^ 

x (l + ~*< 8 >(s)(a fc - x) 3 + i*W(x)(x fc - x) 4 + [ ^ 2 ( 2 )]2 (gfc - ^) 6 



H 1 "i^) +0(1/ ^ ) ;^- (454) 

For the expansion of /(xfe) around X& = x, it is interesting to notice that 

f\x) = 0{l/x 2 ) , and f"(x) = 0{l/x 3 ) . (4.55) 

In the integral, the factor containing f'(x) is multiplied (x k — x) which after integration 
vanishes. Therefore, the relevant contributions to order 0(1/ x 2 ) are 



n hs ( x ) = U l_ [ °° c -ap'(x)(x k -x) 2 /x J °rf( x ) 

X 



X 



(l + $*®(x)(x k - x) 3 + ±¥*>(x)(x k - x) 4 + [ ^ ]2 ( f\ x k - x) 6 ) 



x ( 1 1 —— ) dx k + 0(l/x 3 ) . (4.56) 



\2ap(x 

Then, performing the Gaussian integrals and replacing the dominant values of \l/(x) and its 
derivatives from Eqs. (14. 53p for x ^> 1, we find 

n(x) =n b (l + —^—) (1 - — ^ ) + 0(l/x 3 ) = n fe + 0(l/x 3 ) . (4.57) 



12ax 2 / \ 12ax 2 ^ 

In the bulk of the plasma, the density of particles equal the bulk density, as expected. The 
above calculation, based the Laplace method, generates an expansion in powers of 1/x for 
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the density. The first correction to the background density, in 1/x 2 , has been shown to be 
zero. We conjecture that this is probably true for any subsequent corrections in powers 1/x 
if the expansion is pushed further, because the corrections to the bulk density are probably 
exponentially small, rather than in powers of 1/x, due to the screening effects. In the 
following subsections, we consider the expansion of the density in other types of limits, and 
in particular close to the boundaries, and the results suggest that our conjecture is true. 



8. Asymptotics of the density close to the boundary in the limit x m — > oo 

We study here the density close to the boundary x = x m in the limit x m — > oo and M 
fixed. Since in this limit this region is almost flat, one would expect to recover the result for 



23(. Let x = x m + y where y <C x m is of order 1. 



we have 



the OCP in a flat space near a wall 

Using the dominant term of the asymptotics (14. 15ft . 

BN{k) = ^V^pW) e 2kl ^~ ah ^ [1 + erf (e fe )] , (4.58) 

r, 2k(lnx-lnx k )~a[h(x)-h(x k )] 

n hs (x) = ^±y e (4.59) 

V 71 t^ ^axkp'ixk) [1 + erf (e*)] 

where we recall that Xk = p~ 1 {k/a). The exponential term in the sum has a maximum when 

Xk = x i.e. k = A; max = ap(x), and since x is close to x m — > oo, the function is very peaked 

near this maximum. Thus, we can use Laplace method to compute the sum. Expanding the 

argument of the exponential up to order 2 in k — A; max , we have 

baf \ 2n b ^-4 eX P ~axp'(x)(k fcn — j 

n (x) = > . (4.60) 

V^t^o vWR[l + erf(e fc )] 

Now, replacing the sum by an integral over t — and replacing x = x m — y, we find 

_ 2n b r ex p [-(* - ^y) 2 ] 



Jo 1 + erf(t) 

Since both x m — > oo, and x — > oo, in that region, the space is almost flat. If s is the 



geodesic distance from x to the border, then we have y ~ a/ (nn b /a) s, and equation (I4.6ip 
reproduces the result for the flat space [23(, as expected. 
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9. Density in the thermodynamic limit at fixed shape: a — > oo and x m fixed. 
Using the expansion (1B18I) of B^{k) for the fixed shape situation, we have 

N_Z} „-a[h(x)—2p(x k ) lnx-h(x k )+2p(x h )lnx k ] 

n hs (x)=2n b Y] = . (4.62) 

y/anx k j/(x k ) [erf(e M ) + erf(e fc , m )] 

Once again, to evaluate this sum when a — > oo it is convenient to use Laplace method. The 
argument of the exponential has a maximum when k is such that x k = x. Transforming the 
sum into an integral over x k , and expanding the argument of the integral to order (x k — x) 2 , 
we have 

h A 2n b ^/a f Xm lp'(x k ) e -<*p'(x)(x-*k) 2 /x 



n *( z ) = I d$k (463) 

V 71 " Ji V Xk er *( e M) + erf(e fcim ) 
Depending on the value of x the result will be different, since we have to take special care 
of the different cases when the corresponding dominant values of x k are close to the limits 
of integration or not. 

Let us first consider the case when x — 1 and of order one. This means we are 

interested in the density in the bulk of the system, far away from the boundaries. In this 
case, since e^i and e k>mi defined in (1B19j) and (1B20j) . are proportional to yfot — > oo, then 



each error function in the denominator of (I4.63P converge to 1. Also, the dominant values 
of x k , close to x (more precisely, x — x k of order l/^/a), are far away from 1 and x m (more 
precisely, x k — 1 and x m — x k are of order 1). Then, we can extend the limits of integration 
to — oo and +oo, and approximate x k by x in the term p'(x k )/xk- The resulting Gaussian 
integral is easily performed, to find 

n(x) = rif, , when x — 1 and x m — x are of order 1. (4.64) 

Let us now consider the case when of order 1/y/a, i.e. we study the density 

close to the boundary at x m . In this case e k , m is of order 1 and the term eri(e k ,m) cannot 
be approximated to 1, whereas e^i oc yfa — > oo and erf^^i) — > 1. The terms p'(x k )/x k and 
p'(x)/x can be approximated to p'(x m )/x m up to corrections of order l/y/a. Using t = e k ^ m 
as new variable of integration, we obtain 

2" 



n s {x) = -= \ 
Jo 



cxp 

-oo 



■ ' \ I ~ ^ (Xrn. X j t 

dt , for x m — x of order . 
erf(t) sjot 

(4.65) 




FIG. 2: The normalized one-body density n hs (x)/rib, for the 2dOCP on just one universe of the 
surface S. The dashed line corresponds to a numerical evaluation, obtained from (|4.4ip . with 
N = 100, x m — 2 and a — 4.15493. The full line corresponds to the asymptotic result in the fixed 
shape limit when a — > oo, and x m = 2 fixed. 

In the case where x — 1 is of order 1/a/«, close to the other boundary, a similar calculation 
yields, 

2" 



n hs (x) 



2nb 

^ JO 

where p'(l) = 32. 



exp 



t- v^(T)(x-i))' 



1 + erf(t) 



— dt , for x — 1 of order 



1 



(4.66) 



Fig. [2] compares the density profile for finite iV = 100 with the asymptotic results (14.641) . 
(I4.65P and (14.661) . The figure show how the density tends to the background density, n b , 
far from the boundaries. Near the boundaries it has a peak, eventually decreasing below 
rib when approaching the boundary. In the limit a — > oo, the value of the density at each 
boundary is nbln2. 

Interestingly, the results (14.641) . ( 14.651) and (14.661) turn out to be the same than the one 
for a flat space near a hard wall 23j. From the metric (12.161) . we deduce that the geodesic 



distance to the boundary at x m is s = M(x m - x) ^ p' (x m ) / (8x m ) (when of order 
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1/y/a), and a similar expression for the distance to the boundary at x = 1 replacing x m by 
1. Then, in terms of the geodesic distance s to the border, the results (I4.65P and 



exactly the same as those of an OCP in a flat space close to a plane hard wall 231 ] . 



4.66P are 



n(s) 



2nb 
ft 



exp 



dt. 



(4.67) 



1 + erf (t) 

This result shows that there exists an interesting universality for the density, because, 
although we are considering a limit where curvature effects are important, the density turns 
out to be the same as the one for a flat space. 



B. The 2dOCP on the whole surface with potential — ln(|z — z'\/ yf\zz'\) 

1. Partition function 

Until now we studied the 2dOCP on just one universe. Let us find the thermodynamic 
properties of the 2dOCP on the whole surface S. In this case, we also work in the canonical 
ensemble with a global neutral system. The position Zk = x^ lVk of each particle can be in 
the range l/x m < x^ < x, m . The total number particles N is now expressed in terms of the 
function p as iV = 2ap(x m ) . Similar calculations to the ones of the previous section lead to 
the following expression for the partition function, when (3q 2 = 2, 

now, with 



Z ws exp(-/?F WS ) (4.68) 



N 2 r Xm \n(r)V 

- /3F WS = Nb + Nah{x m ) In x m - a 2 / dx (4.69) 

2 Jl/Xm. X 



1/Xr, 

and 

N 

|2 



z r = ^ [ U dS ^~ ah{x ^ N+1 n 

' "* i=l Ki<i< 



Zi-zj\ 2 . (4.70) 

l<i<j<N 

Expanding the Vandermonde determinant and performing the angular integrals we find 

jV-l 

= n ^ ( 4 - 71 ) 

with 

B N (k) = J x ^-N+l e -ah{x) dS (4>72) 

x 2k ~ N+l e- ah{x) p\x)dx. (4.73) 



a 
n 



l/Xn 
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The function B^{k) is very similar to B^, and its asymptotic behavior for large values of N 
can be obtained by Laplace method as explained in appendix [0 

2. Thermodynamic limit R —* oo, x m — * oo, and fixed M 



Writing the partition function as 



A' 



In Z WS = ln BnW ~ ln ®n{N) 



k=0 



and using the asymptotic expansion ( 1B31I) for Bn, we have 



InZ™ 



ln + + s 2 ws + s™ + sr + sr 



2tt 

ln x m - N ln x m + ah(x m ) 



ln 



where 



A? 



5 wi 



s ws 



s ws 



5> 



fc=0 

A' 



aa^AL 1 + 



l — n V / 



ah(x k 



k=0 

N 

fc=0 
N 

fc=0 

£6n 



erf(e fc>min ) + erf(e fc 

,: 



fc'=l 



12 



8/ A/ 



-i 

E 

k'=-N/2 



1 

12 



MA 

8/ A;' 



6 



lnx m + C>(l) 



(4.74) 



(4.75) 



(4.76) 

(4.77) 
(4.78) 
(4.79) 
(4.80) 



and e^min and efc max are defined in equation flB33j) . Notice that SJ S = due to the symmetry 
relation x_£ = l/xe, therefore only the sums 5f s , S™, SJ S and S™ s contribute to the result. 
These sums are similar to the ones defined for the half surface case, with the difference 
that the running index k' = k — N/2 varies from —N/2 to N/2 instead of to iV as in 
the half surface case. This difference is important when considering the remainder terms in 
the Euler-McLaurin expansion, because now both terms for k' = —N/2 and k' = N/2 are 
important in the thermodynamic limit. In the half surface case only the contribution for 
k = N was important in the thermodynamic limit. 
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The asymptotic expansion of each sum, for x m — > oo, is now 
S™ = ^lna + x^(21nx m -l) + 2x m (81nx m -4) + (28a + l)lnx m + 12a(lnx m ) 2 + 0(l) 



JV 2 P™- [p(x)l 2 1 
= -^-lnx m + a 2 / dx — aNh(x m ) + N\nx m — ah(x m ) + — \nx m + 0(1) 

2 Jl/Xm. x 3 



(4.81) 

'I/Xr, 

(4.82) 

/ 47TQ; 

Sr = -2x m J /5 7har d + 0(l) (4.83) 

V n b 

where 7 is defined in equation (14.281) . The free energy is given by (3F WS = — lnZ ws , with 



In Z ws = 2ax 2 m In x m + N ( b + In ) - ax m + 8ax m (2 In x m - 1) - 2C fi /37 har d 



+12a(lnx m ) 2 + 28alnx m + - lnx m + 0(1) . (4.84) 

6 

We notice that the free energy for this system turns out to be nonextensive with a term 



2x^ n \nx m . This is probably due to the special form of the potential — \n(\z — z'\/^/\zz'\): 
the contribution from the denominator in the logarithm can be written as a one-body term 
{(N — l)/2] lnx, which is not intensive but extensive. However, this nonextensivity of the 
final result is mild, and can be cured by choosing the arbitrary additive constant bo of the 
Coulomb potential as 60 — — hi(Mx m ) + constant. 

3. Thermodynamic limit at fixed shape: a — > 00 and x m fixed 



For this situation, we use the asymptotic behavior (1B34j) of B 



>JV 



In Z WS = N In + SJ sMcd + SJ sfixcd + SJ sfixcd + SJ sMed (4.85) 
rib 



where, now 



JV-l 



sr Med = ^E ln [Vfp'(Vf)i (4-86) 

fc=0 
JV-l 

^ws.fixed = - a J2[h(x k _N) -2p(x k _N)\nx k _N} (4.87) 

fc=0 



gws.fixed _ erf(efc,min) + erf(efc, m ax) ^ gg^ 



k=0 
JV-l 



5 w S) fixed = J2\nx k _N (4.89) 



2 

k=0 
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These sums can be computed as earlier using Euler-McLaurin summation formula. We 
notice that 

S 4 ws ' fixed = a / \nxp'(x)dx + 0(l) = + 0(1) (4.90) 

J 1/ x m 

because of the symmetry properties ln(l/x) = — lnx and p'(l/x)d(l/x) = —p'(x)dx. In the 
computation of 5 , ™ s,fixed there is an important difference with the case of the half surface 



section, due to the contribution when k = 0, since x_n/2 = l/^jv/a — V x r, 

There is no 0(a) contribution from g^ s ^ xed _ Finally, the free energy /3F WS = — lnZ ws is 
given by 



Sr Med = -aNh{x m ) -—\nx m + a 2 H dx + 0(1) . (4.91) 

2 J 1/Xrn x 



\nZ ws = a 



A% b J J 1/Xm X 3 X 2 

2 



/ Ana I I 

-2W X m 1 + — ) /5 7 hard + 0(l) (4.92) 

We notice that the free energy has again a nonextensive term proportional to a In a, but, 
once again, it can be cured by choosing the constant b as b = — ln(Mx m ) + constant. 
The perimeter correction, 2Cij/?7h ar d? proportional to y/a, has the same form as for the half 
surface case, with equal contributions from each boundary at x = l/x m and x = x m . Once 
again, there is no In a correction in agreement with the general theory of Ref. 0, Q and the 
fact that the Euler characteristic of this manifold is % = 0. 

4- Density 

The density is now given by 

N ~ 1 r 2k-N+l p -ah(x) 

- WS ^) = E — s-^ — ( 4 - 93 ) 

Due to the fact that the asymptotic behavior of B^ik) is almost the same as the one of 
B]sj{k') with k! = \k — y|, the behavior of the density turn out to be the same as for the half 
surface case, in the thermodynamic limit a — > oo, x m fixed, 

n(x) = rib , in the bulk, ie., when x — x m and x are of order 1. (4.94) 

Xm. 
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And, close to the boundaries, x — > Xb with Xb = x m or Xb = l/x r , 

21 



n\x) 



2n b 



+ OC 



exp 



-It 



— dt , for Xb — x of order —= . (4.95) 



1 + erf(t) 

If the result is expressed in terms of the geodesic distance s to the border, we recover, once 
again, the result of the OCP in a flat space near a hard wall ( 14. 67ft . 



C. The 2dOCP on the half surface with potential - ln(\z - z'\/y/\zz'\) 

1. Partition function 

In this case, we have N = ap{x m ). Following similar calculations to the ones of the 
previous cases, we find that the partition function, at [3q 2 = 2, is 

= Zfe-^F (4.96) 



with 



and 



with 



/3F hs = a 2 p(x m )h(x m ) - p(x m ) 2 \nx m + \ " M^IL dx _ Nbo (4.97) 

x 



N-l 



fc=0 



n r Xm 

B N (k) = — x 2k+1 e~ ah{x) dx (4.99) 
nb Ji 

2. Thermodynamic limit R — > 00 ; x m — > 00, and fixed M 

The asymptotic expansion of B^{k) is obtained from equation (1B31I) replacing k' by fe 
and considering only the case k > 0. As explained in appendix [Bj the main difference with 
the other half surface case (section II V Al) . is an additional term Xk in each factor of the 
partition function and the additional term (3 /{8k)) in the expansion ( 1B31I) . Therefore, the 
partition function can be obtained from the one of the half surface with potential — In \z — z'\ 
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by adding the terms 

N-l 

Sf= ^lna; fc , (4.100) 

k=0 

5 ? = E^ = ^ lniV + ( 1 ) = i lna; - +0 ( 1 )- ( 4 - 101 ) 

k=l 

Using Euler-McLaurin expansion, we have 

N 

5*4 S = y]ln^ fc - In x m 

/•Xm. J 

= / ap'(x) laxdx + — \nx m — \nx m + 0(1) 
Ji 2 

= ap(x m ) lna; m - a j 



P ^ dx-\\nx m + 0(l) 



x 2 

ap(x m ) In x m - -ah(x m ) - - In x m + O(l) , (4.102) 



where we used the property (13.201) . Finally 



\nZ hs = ax 2 m \nx m + N \ b + In— ) - %-x 2 m + Aax m (2\nx m - 1) 

\ \ 2 n b J 2 

-C R /3 7hard + 6a(lnx m ) 2 + 14alnx m + ^lnx m + 0(1) . (4.103) 
The result is one-half of the one for the full surface, lnZ ws , as it might be expected. 



3. Thermodynamic limit at fixed shape: a —* oo and x m fixed 
For this case, the asymptotics of Bn are very similar to those of Bn from equation flB18j) 



B N (k) ~ x k B N {k) . (4.104) 

Therefore, the only difference from the calculations of the half surface case with potential 
— In \z — z'\ + constant, and this case, is the sum 

_ JV-l 

^hs.fixed = J2\nx k . (4.105) 

fc=0 

We have 



S±' xe = / ap'(x) Inxdx + 0(1) 



ap(x m ) \nx m - \ah(x m ) + 0(1) . (4.106) 
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Here, the term k = N and the remainder of the Euler-McLaurin expansion give corrections 
of order O(a ) = 0(1), as opposed to the previous section where they gave contributions of 
order 0(lnx m ). 



Finally, we find 



a 



p(x r 



1 v / 2mr , , 

2 ln — — + bo} + 



4na 
n b 



1 



n b 



1 + — ) +4 

x, 



°°(l + x)\ (1 + x) 4 J 
^ — in - — dx 



X" 



x- 



/37hard + 0(l). 



(4.107) 



The bulk free energy, proportional to a, plus the nonextensive term proportional a In a, are 
one-half the ones from equation (I4.92p for the full surface expected. The perimeter 

contribution, proportional to \fa is again the same as in all the previous cases of thermody- 
namic limit at fixed shape, i.e. a contribution /37h ar dC& for each boundary at x b = x m and at 
x\, = 1 with perimeter C b (14.381) . Once again, there is no In a correction in agreement with 
the fact that the Euler characteristic of this manifold is \ = 0. 



D. The grounded horizon case 

1. Grand canonical partition function 

In order to find the partition function for the system in the half space, with a metallic 
grounded boundary at x = 1, when the charges interacting through the pair potential of 
Eq. (I3.12p it is convenient to work in the grand canonical ensemble instead, and use the 



techniques developed in Refs. [4j, |25(]. We consider a system with a fixed background density 
p b . The fugacity ( = e^/A 2 , where /i is the chemical potential, controls the average number 
of particles (N), and in general the system is nonneutral (N) ^ N b , where N b = ap(x m ). 
The excess charge is expected to be found near the boundaries at x = 1 and x = x m , while 
in the bulk the system is expected to be locally neutral. In order to avoid the collapse of a 
particle into the metallic boundary, due to its attraction to the image charges, we confine 
the particles to be in a "disk" domain CIr, where x G [1 + w,x m ]. We introduced a small 
gap w between the metallic boundary and the domain containing the particles, the geodesic 
width of this gap is W = ^ap'(l)/ {2-nn b ) w. On the other hand, for simplicity, we consider 
that the fixed background extends up to the metallic boundary. 
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In the potential energy of the system (14. ip we should add the self energy of each particle, 
that is due to the fact that each particle polarizes the metallic boundary, creating an induced 
surface charge density. This self energy is \ ln[|x 2 - 1\M/2L], where the constant ln(M/2L) 
has been added to recover, in the limit M — > 0, the self energy of a charged particle near a 
plane grounded wall in flat space. 

The grand partition function, when (3q 2 = 2, is 



gll 



oo . N N 



N=l 



i=l 



i<3 



Zi Zi 



1 ZiZj 



2 N 

n 

i=i 



i^p-ir 1 



N 

n 

i=l 



-a[h(xi)-2N b \nxi] 



(4.108) 

where for N = 1 the product Yli<j mus t be replaced by 1. The domain of integration for 
each particle is Qr. We have defined a rescaled fugacity ( = 2L(/M and 



- (3Ft = aN b h(x m ) - Nl lnx m - a 2 J 



dx 



(4.109) 



which is very similar to Fq s , except that here N b = ap(x m ) is not equal to N the number of 
particles. 

Let us define a set of reduced complex coordinates Ui = Z{ and its corresponding images 
u* = By using Cauchy identity 



det 



1 



1 1 n w («i-«5) 



(4.110) 



the particle-particle interaction and self energy terms can be cast into the form 

2 N 



n 

Kj 



Z', Zq 



1 Zi Zj 



ZiZ j / (i,i)e{i,-,7V} 2 



(4.111) 



The grand canonical partition function is then 

N N 



c = e 



x,)l det 



1 



I . , i - ( 4 - 112 ) 

N=l " i=l i=l ^ 1 i / _ 

with C(x) = ^e - "^^ -2 ^ 111 ^. We shall now recall how this expression can be reduced to a 
Fredholm determinant j^]. Let us consider the Gaussian partition function 



VipVip exp 



(4.113) 



The fields ■?/> and ip are anticommuting Grassmann variables. The Gaussian measure 
in (14.1131) is chosen such that its covariance is equal to 

1 



(^(q i )^(q J -)> = A(zi,Zj 



1 Zi Zj 



(4.114) 
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where (...) denotes an average taken with the Gaussian weight of (14.1131) . By construction 
we have 

Z = detiA- 1 ) (4.115) 
Let us now consider the following partition function 



VipVip exp 
which is equal to 



ij(<\)A-\z,z')^{c{)dSds' - / c(aOV<qMq) dS 



(4.116) 



Z = det(A- 1 -() (4.117) 

and then 

— = det[A(A- 1 - C)] = det(l + K) (4.118) 
where K is an integral operator (with integration measure dS) with kernel 

K(q, q') = -C(x') A(z, z') = ■ (4.H9) 

Expanding the ratio Z/Z in powers of ( we have 

„ oo N N 

Y = 1 + T, m / II ^(- 1 ) 7V II CO*) <^(qi)^(qi) • • ■ ^(q* Mq*)> (4.120) 

jv=1 ' J i=l i=l 



Now, using Wick theorem for anticommuting variables (26j, we find that 

(^(qi)V(qi)'-"0(qjv)^(qiv)) = detA(z i ,^) = det ( 1 - ) (4.121) 

\1 z i z j / 

Comparing equations (14.1201) and (14.1 12[) with the help of equation (14.1211) we conclude that 

H = e -f>*t JL = e -0*f det (l + X) (4.122) 

The problem of computing the grand canonical partition function has been reduced to 
finding the eigenvalues A of the operator K. The eigenvalue problem for K reads 

-/ ^-^®{x',y')dS' = A$(x,v?) (4.123) 
Jn R 1 ~~ zz 

For A 7^ we notice from equation (14.1231) that <&(x,<p) = *&(z) is an analytical function 
of z — xe ltp in the region \z\ > 1. Because of the circular symmetry, it is natural to try 
<&(z) = $i(z) = z~ l with I > 1 a positive integer. Expanding 

j oo 



1 - 

n=l 



(4-124) 
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and replacing $e(z) = z in equation (14.1231) . we show that is indeed an eigenfunction 
of K with eigenvalue 

\ i = {B$ i (N b -e) (4.125) 

where 

Bf b (k) = — x 2k e- ah{x) p'{x)dx (4.126) 

which is very similar to Bn defined in Eq. (14. lip , except for the small gap w in the lower 
limit of integration. So, we arrive to the result for the grand potential 



pn 



oo 

lnS = /3F| h -^ln [l + (Bf b (N b 



(4.127) 



2. Thermodynamic limit at fixed shape: a — > oo and x m fixed 



Let us define k = N b — i for i e N*, thus k is positive, then negative when i increases. 

Therefore, it is convenient to split the sum (14.1271) in In H into two parts 

-i 

£ ln[l + CB? 6 (*)] (4-128) 



fixed 



fc=— oo 
AT 6 -1 



(4.129) 



fc=0 



The asymptotic behavior of 13$ (k) when a — > oo can be directly deduced from the one of 
Sat found in appendix [Bj Eq. (1B18j) . taking into account the small gap w near the boundary 
at x — 1 + w. When k < 0, we have x^ < 1, then we notice that e^i defined in (IB20p 
is negative, and that the relevant contributions to the sum Sf 11 ' ed are obtained when k is 
close to 0, more precisely k of order 0(y/Nb). So, we expand Xk around Xk = 1 up to order 
(xk — l) 2 in the exponential term e~ a ^- h( - Xk ^~ 2p( - Xk ^ lnx ^ from Eq. flB18j) . Then, we have, for 
k < of order 0(y/N~ b ), 



Bf b (k) = VWW e ^(D fl-* fe ) 2 erfct^l) (1 + w - x h )) 
z,n b 



(4.130) 



where erfc(w) = 1 — erf(tt) is the complementary error function. Then, up to corrections 
of order 0(1), the sum 5'| h ' fixed can b e transformed into an integral over the variable t = 
\J ap'(l) (1 — Xk), to find 



Si 



;h, fixed 



'ap'{l) / In 
'o 



2n h 



-e* erfc (t + V2lF^ b W) 



eft + 0(1). (4.131) 
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Let C\ = ^/27rap'(l)/n6, be total length of the boundary at x — 1. We notice that 

y/a-Kp'{\) _ (d 2(L d 



c- 



2n b 



/ 2n h 



/ 2~n h M 



(4.132) 



is fixed and of order 0(1) in the limit M — > oo, since in the fixed shape limit C\/M is fixed. 
Therefore S$' 6 gives a contribution proportional to the perimeter C\. 



For Sf ' xed , we define 



and we write 

^gh, fixed 



i = \/ cnp'(l) (1 + W - X k ) 



(4.133) 



1 + C^/airxkp'ixk) c _ a[fe ( £fc )_ 2p(£fc ) ln£fc1 



fc=0 

r<sh, fixed i rihs, fixed . nhs, fixed . tvt i <^^/d7T 

51 +5, + 6 2 +A/ 6 ln 

n 6 



erf(e M ) + erf(e fc . 



(4.134) 



where 



St 



fixed 



N b -l 
fc=0 



3 a[/i(x fc )-2p(£ fc )ln£ fc ] 



(4.135) 



and we see that the sums S± s ' xed and g^' 6 * 6 ' 1 reappear. These are defined in equations (I4.3ip 
and (14.321) and computed in f)4.34p and (|4.35p . In a similar way to 5 r | h ' flxed j ^-^fixed gj ves on iy 
boundary contributions when k is close to 0, of order (grounded boundary at x — 1) 
and when k is close to ]V& with N b — k of order V^fe (boundary at x = x m ). We have, 



fixed 



VV(1) / In 



n 6 e 



/ In 

'o 



erf(t) + 1 



tit 



+ [erf( t _ ^T bW ) + 1] 
'(1) 1 J 

dt (4.136) 



Let us introduce again the perimeter of the outer boundary at x 
2-KaXmp' (x m ) / rib. Putting together all terms, we finally have 



In! 



-N b f3u B + - [h(x m ) - 2p(x m ) \nx m ] + a — In — dx 

2 J 1 x 6 x 2 

-Ci/37 meta i - C iJ /?7 hard + 0(1) 



(4.137) 



where 



f3u)B = — In 



2nCL 
\/2ru 



(4.138) 
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is the bulk grand potential per particle of the OCP near a plane metallic wall in the flat space. 
The surface (perimeter) tensions 7 me tai and 7h ar d associated to each boundary (metallic at 
x b — 1, and hard wall at Xb = x m ) are given by 



metal 





In 



In 



2n b 



dt 



n b e 



-f 



. r -= + ~ [erf (t - V2^T b W) + 1] 



dt (4.139) 



with x b = 1, and (I4.28P for /?7hard- 

Notice, once again, that the combination 



(y/airxhp'ixh) _ 2(L C b 



(4.140) 



2n b V2rT b M 

is finite in this fixed shape limit, since the perimeter C b of the boundary at x b scales as M. 
Up to a rescaling of the fugacity ( to absorb the factor C b /M, the surface tension near the 
metallic boundary 7 me tai is the same as the one found in Ref. 
similar to the one found in Ref. 



4| in flat space. It is also 



25] with a small difference due to the fact that in that 



reference the background does not extend up to the metallic boundary, but has also a small 
gap near the boundary. 

There is no In a correction in the grand potential in agreement with the fact that the 
Euler characteristic of the manifold is \ = 0. 

Let us decompose In S into its bulk and perimeter parts, 



InS = -PQf - Ci/37 m etal - Cif^hard + 0(1) 



(4.141) 



with the bulk grand potential Qf given by 



(3fLf = -N b (3ujB + ^ [h(x m ) - 2p(x m ) hxx m \ + a 



a 



[l + x)\ (1+x) 4 , 
^ln^ ^dx. (4.142) 



x" 



ar 



The average number of particles is given by the usual thermodynamic relation (N) = 
(d(\nE)/d(. Following (14.1411) . it can be decomposed into bulk and perimeter contribu- 
tions, 

(N) = N b - CiC^lp ■ (4.143) 

The boundary at x = x m does not contribute because 7hard does not depend on the fugacity. 
From this equation, we can deduce the perimeter linear charge density a which accumulates 
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near the metallic boundary 

a = -C ^ cta ' • (4.144) 

We can also notice that the bulk Helmoltz free energy Ff 1 = flf 1 + /xAT b is the same as for 
the half surface, with Coulomb potential G hs , given in ( 14.371) . 



3. Thermodynamic limit R — > oo, x m — > oo, and fixed M 



This limit is of restricted interest, since the metallic boundary perimeter remains of order 
0(1), we expect to find the same thermodynamic quantities as in the half surface case 
with hard wall "horizon" boundary up to order 0(hxx m ). This is indeed the case: let us 
split InH into two sums Sq and Sf 1 as in ( I4.128P and (14.1291) . For k < 0, the asymptotic 
expansion of Bx b (k) derived in appendix [B] should be revised, because the absolute maximum 
of the integrand is obtained for values of the variable of integration outside the domain of 
integration. Within the domain of integration the maximum value of the integrand in (14. 1 26[) 
is obtained when x = 1 + w. Expanding the integrand around that value, we obtain to first 
order, for large \k\, 



2n b \k\ 



(4.145) 



Then 



D 6 



ln[l + 0B&(*) 

fc=— oo 

POO 

dkln 

'o 

0(1), 



ap'(l + w) 2wM 
^ 2n b \k\ 



+ 0(1) 



(4.146) 



does not contribute to the result at orders greater than 0(1). For the other sum, we have 



k=0 k=0 

= ^ In \(Bf b (k)} +0(1). 



C<(fc) 



k=0 



(4.147) 



The second sum is indeed 0(1), because l/[(B^ b (k)] has a fast exponential decay for large 
k, therefore the sum can be converted into an finite [order 0(1)] integral over the variable 
k. 
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Now, since the asymptotic behavior of B^(k), for k > and large, is essentially the same 
as the one for £>jv 6 (/c), we immediately find, up to 0(1) corrections, 

InH = (3fiN b + \nZ hs + 0(1) (4.148) 

where In Z hs is minus the free energy in the half surface case with hard wall boundary, given 
by gZSD. 



4- The one-body density 

As usual one can compute the density by doing a functional derivative of the grand 
potential with respect to a position-dependent fugacity C(q) 

™ gh (q) = C(q)^- (4-149) 

For the present case of a curved space, we shall understand the functional derivative with 
the rule = <5(q, q') where 5(q, q') = 5(x — x')5(<p — ip')/^/g is the Dirac distribution on 
the curved surface. 



Using a Dirac-like notation, one can formally write 

In 5 = Tr ln(l + K) - /5F gh = J (q |ln(l - C{q)A) | q) dS - /3F gh 
Then, doing the functional derivative (14. 149ft . one obtains 



(4.150) 



rr 



q) = C <q I (1 + K)-\-A) | q> = CG(q, q) (4.151) 



where we have defined G(q, q') by G = (1 + K) 1 {— A). More explicitly, G is the solution 
of (1 + K)G = -A, that is 

0(q, q') - / C(x")^$- dS" = . (4.152) 

' ' 1 — zz" 1 — zz' 



R 



From this integral equation, one can see that G(q, q') is an analytical function of z in the 
region \z\ > 1. Then, we look for a solution in the form of a Laurent series 

oc 

G(q,q') = J>^ r V* (4-153) 

into equation (14.1521) yields 
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Recalling that X £ = (B^(N b — £), the density is given by 

^b -1 2k -ah(x) 

n*(x) =(V (4.155) 

5. Density in the thermodynamic limit at fixed shape a — > oo and x m fixed. 

Using the asymptotic behavior ( 1B18I) of Bff , we have 

n gh^ - ( ^ ex P (- a lH x ) ~ 2p(x fc ) \nx - h(x k ) + 2p(x k ) \nx k }) ^ 

fe= _ oc e alh(x h )-2 P (x k )lnx k ] + Sy/™xrf{x k ) foffa^ + e rf(e fc , m )] 

Once again, this sum can be evaluated using Laplace method. The exponential in the 
numerator presents a peaked maximum for k such that x k = x. Expanding the argument of 
the exponential around its maximum, we have 

N b -ap'(x)(x-x k ) 2 /x 

n gh (x) = C Yl 7T=W^ • ( 4 - 157 ) 

e a[h(x k )-2 P (x k )inx k ] + ^V»»^ [ er f(e M ) + erf(e fc>ro )] 

Now, three cases has to be considered, depending on the value of x. 

If x is in the bulk, i.e. x — 1 and a; m — a; of order 1, the exponential term in denominator 
vanishes in the limit a — > oo, and we end up with an expression which is essentially the same 
as in the canonical case f)4.62p [the difference in the lower limit of summation is irrelevant 
in this case since the summand vanishes very fast when x k differs from x}. Therefore, in the 
bulk, n gh (x) = n b as expected. 

When of order 0(l/y/a), once again the exponential term in the denominator 

vanishes in the limit a — > oo. The resulting expression is transformed into an integral over 
the variable e^m, and following identical calculations as the ones from subsection IIV A 9l we 
find that, n gh (x) = n hs (x), that is the same result (14.651) as for the hard wall boundary. This 
is somehow expected since, the boundary at X — OC is of the hard wall type. Notice that 
the density profile near this boundary does not depend on the fugacity £. 

The last case is for the density profile close to the metallic boundary, when x — 1 is of 
order 0(l/^/a). In this case, contrary to the previous ones, the exponential term in the 
denominator does not vanish. Expanding it around x k = 1, we have 

N b -ap'(x)(x-x k ) 2 /x 

n gh (x) =( V . (4.158) 

^oo e -<i + iv^gM [erf( g M ) + i] 
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Transforming the summation into an integral over the variable t = — 1, we find 



,"+°° -[t+y/ap'(l)(x-l)] 2 J, 

n*\x) = Cv^O) / 7^== • (4-159) 

> e -i 2 + ^erfc(t + ^) 



For purposes of comparison with Ref. 



25], this can be rewritten as 



2y/ aP '(l)(x-l)t ^ 



/(x) = ( v / «p ; (l)e- y(1)[(l - 1 -'" ) -'" 1 j- . (4.160) 

■/-«> x + Vgw erfc ( f ) e (t-V2^ty) 2 

Which is very similar to the density profile near a plane metallic wall in flat space found 



in Ref. 25J [there is a small difference, due to the fact that in 25J the background did 
not extend up to the metallic wall, but also had a gap, contrary to our present model]. 
Fig. [3] shows the density profile for two different values of the fugacity, and compares the 
asymptotic results with a direct numerical evaluation of the density. 

Interestingly, one again, the density profile shows a universality feature, in the sense that 
it is essentially the same as for a flat space. As in the flat space, the fugacity controls the 
excess charge which accumulates near the metallic wall. Only the density profile close to 
the metallic wall depends on the fugacity. In the bulk, the density is constant, equal to the 
background density. Close to the other boundary (the hard wall one), the density profile 
is the same as in the other models from previous sections, and it does not depend on the 
fugacity. 

V. CONCLUSIONS 

The two-dimensional one-component classical plasma has been studied on Flamm's 
paraboloid (the Riemannian surface obtained from the spatial part of the Schwarzchild met- 
ric). The one-component classical plasma had long been used as the simplest microscopic 



model to describe many Coulomb fluids such as electrolytes, plasmas, molten salts 27]. 
Recently it has also been studied on curved surfaces as the cylinder, the sphere, and the 
pseudosphere. From this point of view, this work presents new results as it describes the 
properties of the plasma on a surface that had never been considered before in this context. 

The Coulomb potential on this surface has been carefully determined. When we limit 
ourselves to study only the upper or lower half parts (S±) of the surface (see Fig. [[]) the 



] 




FIG. 3: The normalized one-body density n gh (x)/nb, in the grounded horizon case. The dashed 
lines correspond to a numerical evaluation, obtained from (|4.155p . with N = 100, x m = 2 and 
a = 4.15493 and truncating the sum to 301 terms (the lower value of k is —200). The gap close 
to the metallic boundary has been chosen equal to w = 0.01. The full lines correspond to the 
asymptotic result in the fixed shape limit when a —* oo, and x m = 2 fixed. The two upper curves 



correspond to a fugacity given by C,y/a / '(2n&) = QL^J ir/rib = 1, while the two lower ones correspond 
to (L^ir/rib = 0.1. Notice how the value of the fugacity only affects the density profile close to 
the metallic boundary x = 1. 

Coulomb potential is G hs (q, q') = — In \z — z'\ + constant, with the appropriate set of coor- 
dinates (x,tp) denned in section Hi 21 and z = xe tip . When charges from the upper part are 
allowed to interact with particles from the lower part then the Coulomb potential turns out 



to be G ws (q, q') = —\n(\z — z'\/ y/\zz'\) + constant. When the charges live in the upper part 
with the horizon grounded, the Coulomb potential can be determined using the method of 
images form electrostatics. 

Since the Coulomb potential takes a form similar to the one of a flat space, this allows to 
use the usual techniques [2l, |3] to compute the thermodynamic properties when the coupling 
constant Y = Bq 2 = 2. 
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Two different thermodynamic limits have been considered: the one where the radius R of 
the "disk" confining the plasma is allowed to become very big while keeping the surface hole 
radius M constant, and the one where both R — > oo and M — > oo with the ratio R/M kept 
constant (fixed shape limit). In both limits we computed the free energy up to corrections 
of order 0(1). 

The plasma on half surface is found to be thermodynamically stable, in both types of 
thermodynamic limit, upon choosing the arbitrary additive constant in the Coulomb po- 
tential equal to — InM + constant. The system on the full surface is found to be stable 
upon choosing the constant in the Coulomb potential equal to — ln(Mx m ) + constant where 
Xm = (^R+ Vi?-2M)7(2M). 

In the limit R — > oo while keeping M fixed, most of the surface available to the particles 
is almost flat, therefore the bulk free energy is the same as in flat space, but corrections 
from the flat case, due to the curvature effects, appear in the terms proportional to R and 
the terms proportional to In R. These corrections are different for each case (half or whole 
surface) . 

The asymptotic expansion at fixed shape (a — > oo) presents a different value for the bulk 
free energy than in the flat space, due to the curvature corrections. On the other hand, the 
perimeter corrections to the free energy turn out to be the same as for a flat space. This 
expansion of the free energy does not exhibit the logarithmic correction, In a, in agreement 
with the fact that the Euler characteristic of this surface vanishes. 

For completeness, we also studied the system on half surface letting the particles interact 
through the Coulomb potential G ws . In this mixed case the result for the free energy is 
simply one-half the one found for the system on the full surface. 

In the case where the "horizon" is grounded (metallic boundary), the system is studied in 
the grand canonical ensemble. The limit R — > oo with M fixed, reproduces the same results 
as the case of the half surface with potential G hs up to 0(1) corrections, because the effects 
of the size of the metallic boundary remain 0(1). More interesting is the thermodynamic 
limit at fixed shape, where we find that the bulk thermodynamics are the same as for the half 
surface with potential G hs , but a perimeter correction associated to the metallic boundary 
appears. This turns out to be the same as for a flat space. This perimeter correction 
("surface" tension) /37 me tai depends on the value of the fugacity. In the grand canonical 
formalism, the system can be nonneutral, in the bulk the system is locally neutral, and 
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the excess charge is found near the metallic boundary. In contrast, the outer hard wall 
boundary (at x = x m ), exhibits the same density profile as in the other cases, independent 
of the value of the fugacity. This reflects in a perimeter contribution /37hard equal to the one 
of the previous cases. 

The plasma on Flamm' s paraboloid is not homogeneous due to the fact that the curvature 
of the surface is not constant. When the horizon shrinks to a point the upper half surface 
reduces to a plane and one recovers the well known result valid for the one component plasma 
on the plane. In the same limit the whole surface reduces to two flat planes connected by a 
hole at the origin. 

We carefully studied the one body density for several different situations: plasma on half 
surface with potential G hs and G ws , plasma on the whole surface with potential G ws , and 
plasma on half surface with the horizon grounded. When only one-half of the surface is 
occupied by the plasma, if we use G hs as the Coulomb potential, the density shows a peak 
in the neighborhoods of each boundary, tends to a finite value at the boundary and to the 
background density far from it, in the bulk. If we use G ws , instead, the qualitative behavior 
of the density remains the same. In the thermodynamic limit at fixed shape, we find that 
the density profile is the same as in flat space near a hard wall, regardless of the Coulomb 
potential used. 

In the grounded horizon case the density reaches the background density far from the 
boundaries. In this case, the fugacity and the background density control the density profile 
close to the metallic boundary (horizon). In the bulk and close to the outer hard wall 
boundary, the density profile is independent of the fugacity. In the thermodynamic limit at 
fixed shape, the density profile is the same as for a flat space. 

Internal and external screening sum rules have been briefly discussed. Nevertheless, we 
think that systems with non constant curvature should deserve a revisiting of all the common 
sum rules for charged fluids. 
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APPENDIX A: GREEN FUNCTION OF LAPLACE EQUATION 

In this appendix, we illustrate the calculation of the Green function using the original 
system of coordinates (r, ip). The Coulomb potential generated at q = (r, if) by a unit charge 
placed at q = (r , y?o) with r > 2M satisfies the Poisson equation 

AG(r, cp; r , (p ) = -2n5(r - r )5(ip - ip )/ y/g , (Al) 

where g = det(g^) = r 2 / (1 — 2M/r). To solve this equation, we expand the Green function 
G and the second delta distribution in a Fourier series as follows 



oo 



G(r, V ;r , Vo ) = £ e m ^-^g n (r,r ) , (A2) 

n=— oo 

oo 

% - <po) = ^ E ein{ ^ 0) ' ( A3 ) 



n=— oo 



to obtain an ordinary differential equation for g n 
2M\ d 2 (\ M\ d n 2 



dr 2 ^ I r r 2 / dr r 2 



g n (r,r ) = -5(r -r )/y/g . (A4) 



To solve this equation we first solve the homogeneous one for r < r : g n ,-( r , r o) an d r > r : 
g n ,+(j', tq). The solution is, for n^O, 



g n Ar,r )=A n , ± (^ + V^^) 2n + B n ^+V^2My 2n , (A5) 

and, for n = 0, one finds 

9o,± (r, r ) = A ,± + B 0t± ln(>/F + Vr - 2M) . (A6) 

The form of the solution immediately suggest that it is more convenient to work with the 
variable x = Vr - 2M) 2 /(2M). For this reason, we introduced this new system of 

coordinates (x, <p) which is used in the main text. 
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APPENDIX B: ASYMPTOTIC EXPANSIONS OF B N (k), B N (k) AND B N (k) 



1. Asymptotic expansion of B^{k) 



a. Limit N — > oo, x m — > oo, and fixed a 



(Bl) 



Doing the change of variable s = ap(x) in the integral (14. lip , we have 

i r N 

B N {k) = — / x 2k e- ah{x) ds 
n b Jo 

where x is related to the variable of integration s by s = ap(x). The limit k — > oo and 
N ^ oo can be obtained using Laplace method 28(]. To this end, let us write B^ik) as 

jfc f N/k 



B N (k) = - 



n 



where we made the change of variable t = s/k and we defined 



a 



where 



The derivative of 6k is 



4>kif) = 2 lnx — — h(x) 



x = p (kt/a) 



1 dx ^ 

x dt 
2k 



l-t) 
1-t) 



axp'(x) 

where we have used the definition ( 1B4I) of x and the properties (I3.20p of h and p. 
The maximum of 0&(t) is obtained when t — 1. At this point we have 

2k 



*2(1) 



. i + o(l/^) 

ax k p{x k ) 

4k 2 p'(x k ) +x k p"(x k ) 



x 2 k p'(x k ) 3 



2 + 0{l/Vk) 



6k 3 d 
a 3 p'(x k ) dx 



p'(x) + xp"(x) 
x 2 p'(x) 3 



-6 + 0(l/Vk) 



where 



x k =p 1 (k/a) . 



(B2) 



(B3) 
(B4) 

(B5) 
(B6) 



(B7a) 
(B7b) 
(B7c) 

(B8) 
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Expanding (/>k(t) up to order (t — l) 4 , and defining v = y/k\^(T)\ (t — 1), we have 
B N (k) 



-v 2 /2 



X 



<M 3) (i) ^Vi 4) (i) ^ 6 [0l 3) (i)] 



+ 



Let us define 



3!v^|0^(l)| 3 / 2 4!A;|0' fc '(l)l 2 3! 2 2A^(1) 
iV-fc JV-fc 



62(1) 



+ 0(1/VN) 



dv . (B9) 



(BIO) 



X 2A; V2iV 

which is an order one parameter, since we are interested in an expansion for N and k large 
with iV — k of order y/N. Using the integrals 



-*' 2 dv 



erf ( 



e 



/ 

J — c 

f e-^Vtfo = -(2 + e 2 )e- £2 / 2 

•/ — oo 

/ 

J — c 

/ 

7 — c 



e-^ 2 v A dv = 3 



e-*' 2 v 6 dv = 15 



1 + erf | 
1 + erf 



Vv/2 

(A) 



(Bll) 
(B12) 

- e" £2 / 2 e(3 + e 2 ) (B13) 
- e" e2/2 e(15 + 5e 2 + e 4 ) (B14) 



where erf(z) = (2/ y^) J* e " 2 rf-u is the error function, we find in the limit iV — > oo, fc — > oo, 
and finite e^, 



B N (k) 



7T 



[l + erf(e fc )] 



(B15) 



2|02(1)| n 

The functions £i( e fc) an d ^(efe) contain terms proportional e -6 *, from the Gaussian integrals 
above. However, as explained in the main text, these do not contribute to the final result for 
the partition function up to order 0(1), because the exponential term e~ e ^ make convergent 
and finite the integrals of these functions that appear in the calculations, giving terms of 
order 0(1) and 0(1/ y/N) respectively. 



b. Limit N — > oo ; ot — ► oo, fixed x m 

For the determination of the thermodynamic limit at fixed shape, we also need the asymp- 
totic behavior of BN(k) when a — > oo at fixed x m . We write £>tv(^) as 

ry C Xm 

B N (k) = — e- aMx) ~ 2p{£k)lnx] p'(x)dx, (B16) 
nb Ji 
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where we have defined once again x k by k = ap(x k ). We apply Laplace method for a — > oo. 
Let 

F(x) = h(x) -2p(x k )\nx. (B17) 

F has a minimum for x = x k with F"(x k ) = 2p'(x k )/x k . Expanding to the order (x — x k ) 2 the 
argument of the exponential and following calculations similar to the ones of the previous 
section, we find 

B N (k) = V^(^) e - tt[ ^)-2p(*>)in^] [erf(e M ) +erf(e fc , m )] 

x ( 1 + -S (x k ) + 4= Kl,m(efc,m) + Ci.iKi)]) (B18) 

where 



Cfc,™ = \l^-^-(x m -x k ) , (B19) 



e M = y/ap>(l)(x k - 1) . (B20) 

The terms with the error functions come from incomplete Gaussian integral and take into 
account the contribution of values of k such that x m — x k (or x k — 1) is of order l/y/a, or 
equivalently N — k (or k) of order \fN . 

The functions £o(x k ), £i,i( e fc,i)) an< ^ £i,m( e fc,m) can be computed explicitly, pushing the 
expansion one order further. These next order corrections are different than in the previous 
section, in particular (l/a)£ (x k ) ^ I /(12k). 

However, these next order terms are not needed in the computation of the partition 
function at order 0(1), since they give contributions of order 0(1). Note in particular that 
the term £,o(x k )/a gives contributions of order 0(1), contrary to the previous limit studied 
earlier where it gave contributions of order IniV. Indeed, in the logarithm of the partition 
function, this term gives a contribution 

y CoM = 1 f Xm ap '^ x )^ x ) dx + (i) = 0(1) . (B21) 

2. Asymptotic expansions of B^ik) and B^(k) 

To study B N (k), it is convenient to define k' = k — y, then 

n f Xm 

B N (k) = — x 2k 'e- ah(x) xp\x)dx, (B22) 

U b Jl/Xm 
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which is very similar to 

B N (k) = — x 2k e~ ah{x) xp'(x)dx. (B23) 
nb Ji 

changing k' by k, and taking into account the extended domain of integration [l/x m , 1] for 
Bn- As in the previous section, the asymptotic expansions for Bat (A;) and S^ik) can be 
obtained using Laplace method. Notice that for B N (k), k' is in the range [— — , When 
k' < 0, the maximum of the integrand is in the region [l/x m , 1], and when k' > 0, the 
maximum is in the region [1, x m ]. Due to the fact that the contribution to the integral from 
the region [l/x m , 1] is negligible when k! > 0, the asymptotics for B^{k) will be the same 
as those for Biy{k), for k' > 0, doing the change k — > k'. Therefore, we present only the 
derivation of the asymptotics of Bn- 



a. Limit N —* oo ; x m — > do, and fixed a 
We proceed as for B^(k), defining the variable of integration t = ap(x)/k', then 

B N (k) = M xe k '^'®dt (B24) 

n b J n 



2M 



where <fik'(f) is the same function defined in equation (1B3[) . Now we apply Laplace method 
to compute this integral. The main difference with the calculations done for Bn are the 
following. First, taking into account that k' can be positive or negative, we should note that 



-1 + 0(1/W\k'\) k' >0 
<j>Ul) = { WV ' (B25) 

1 + 0(1/^/1^1) k! <0 



2 + 0(1/71^1) k! >0 
-2 + 0(l/v/fF|) k' <0 



^ 3) (1) = { ' (B26) 



<,(4) 



-6 + 0(1/71^1) A;' > 



^(1) = (B27) 



6 + 0(1/71^1) fc' <0 
Second, we also need to expand x close to the maximum which is obtained for t = 1, 

x = x k \l + a(t - 1) + 6(t - l) 2 + 0((t - l) 3 )] (B28) 
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with 



p{xk> 



and 



x k ,p'(x k > 



p(x kl ) 2 p"(x k i\ 
2x fc /p'(x fc /) 3 



l + 0{l/y/\P\) k'>0 
-\ + 0{l/^\F\) k'<0 

-\ + 0{l/^W\) k' >0 
§ + 0(l/vW|) k'<0 



(B29) 



(B30) 



Notice in particular that for the term b, the difference between positive and negative values 
of k! is not only a change of sign. This is to be expected since the function x is not invariant 
under the change x — > 1/x. 

Following very similar calculations to the ones done for Bn with the appropriate changes 
mentioned above, we finally find 



Bnik) = ^-\/ixax k >v'{x k i)e 



1 + 






(A- 











with 



x [erf (e fc>min ) + erf (e fc)inax )] 

k' > 
■| k' < 



+ 



(B31) 



(B32) 



and 



l ap'(x m ) 



iap'(l/x m ) 



\jx T 



x h n - 



(B33a) 
(B33b) 



The dots in ( 1B31I) represent contributions of lower order and of functions of efc jm i n and ejt jmax 
that give 0(1) contributions to the partition function. Comparing to the asymptotics of 13 n 
we notice two differences: the factor x k > multiplying all the expressions and the correction 
c/\k'\. 



b. Limit N — > oo, a — * oo, and fixed x m 

The asymptotic expansion of in this fixed shape situation is simpler, since we do not 
need the terms of order 1/a. Doing similar calculations as the ones done for Bn taking into 
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account the additional factor x in the integral we find 

B N (k) = ^y/^EM e -aM^)-MM^] [er f( efc , mm ) + erf(e fc , max )] . (B34) 
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